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" On some Definite Integrals and a New Method of reducing a 
Function of Spherical Co-ordinates to a Series of Spherical 
Harmonics." By Arthur Schuster, F.B.S. Eeceived May 
30,— Eead June 5, 1902. 

(Abstract.) 

The expansion of a function f(0) of an angle varying between 
and 7T in terms of a series proceeding by the sines of the multiples of 
depends on the fundamental theorem, 

7T 

I* 

I sin ])Q sin qO clO = 0, 

o 

where p and q are integer numbers not equal to each other. Simi- 
larly if P n denotes the zonal harmonic of degree n, p, = cos 0, and 

dp? 

the expansion of a function of in terms of a series of the functions 
Qj depends on the corresponding theorem, 

+i 

Q W*/* = o, 

where i and n are two integer numbers not equal to each other. In 
many practical applications a continuous function is given by means 
of its numerical values at certain points, e.g., for equidistant values 
of 0. 

Such cases present no difficulty when Fourier's analysis is to be 
employed, because there is in that case a summation theorem exactly 
corresponding to the above integration theorem. If be replaced by 
pirjih, where p takes successively the values 1, 2, 3 ... , the equa- 
tion 

p = n—l 

3 sin (ppirjn) sin {pqirjn) = 
P = o 

will hold true. This allows us to determine the coefficients in the 
case of problems in which discontinuous values of the function at 
equidistant points are known (e.g., hourly readings of temperature or 
barometric pressure). If we assume that all Fourier coefficients beyond 
the nth vanish, n equations are obtained, each of which only contains 
one of the unknown quantities. 

If it is desired to expand a function in terms of cosines, a slight 
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modification must be introduced, the summation theorem in that case 
being 

p — n—l 

h + -J eospr cos (pr + 2 cos (ppir/n) cos (pqirjn) = 0, 

p = i 

the first and second terms representing half the value of the product 
for /o = and /o = % respectively. 

There is no corresponding summation theorem in the case of the 
functions Qj , and the application of the method of least squares leads 
to a series of normal equations, each of which contains all the other 
coefficients. This has been one of the great practical difficulties in 
obtaining an expression for the series of spherical harmonics for the 
earth's magnetic potential. 

F. E. Neumann has tried to overcome the difficulty by calculating 
coefficients a v a% • • . % in such a way that 

" S % QZ (/x p ) Qf (/x p ) = 0. 

Here /ii, /x 2 . . . //, g are the quantities for which the values of the 
function to be represented are known. Neumann's process is equiva- 
lent to attaching weights proportional to a p to the different obser- 
vations, a proceeding against which theoretical objections might be 
urged. 

2. The expansion in terms of a series of cosines and sines being so 
much easier than the direct expansion in terms of a series of the 
functions Q£, I have endeavoured to obtain the latter series by means 
of the former. 

It is well known that a function of an angle 0, which is confined 
to the values lying between and ir, may be put either into the form 

a + cii cos 9 + a 2 cos 20 + . . . a p cos p6 + , . . . , 
or into the form 

hi sin + h 2 sin 20 + . . . b p sin p6 + .. . . 

The reduction to the series of spherical harmonics is accomplished 
by calculating and tabulating the coefficients in the series 

gob p0 = A*Q* + A«" +1 QJ +1 + . . . A£Q£ + ...., 

sinpo = iv a q: + B* +1 QJ +1 + . . . B'QJ + . . . 

The choice between the cosine and sine series is open to us, but it 
appears that great simplicity is gained by taking the former series 
when cr is odd and the latter when o- is even. For in that case the 
coefficients A,fand Bj will all vanish, as long as '■ n is smaller than/;. 
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When it is therefore desired to retain only terms as far as the nth. 
degree, the Fourier coefficients need only be calculated as far as 
p = n+1* The position of the earth's magnetic axis, e.g., only 
depending on the terms of the first degree, is completely determined 
by the coefficients b. 2 for a- — and a , a% for or = 1. 

3. The symbolical representation of the results of this paper is much 
facilitated by the introduction of a separate symbol for the product of 
alternate factors, n . n-2 . n-4c . . . 1, if n be odd, or n.n — 2 
. . . 2 if % be odd. I propose to write n 1 ! for such products, and if a 
name be required for the product to call it the " alternate factorial " 
or the " double factorial." Full advantage of the new symbol is only 
gained by extending its meaning to negative values of n. Its complete 
definition may then be included in the equations 

nil = n(n~2)\\, 111-1, 2 !! = 2 . 
From this we may derive when n is negative and odd 

M.4-1 -j 

while for n negative and even, the factorial becomes infinitely large. 

A. The calculation of the factors A % and B.J depends on the values 
of the definite integrals 

-l +i 

Q% cos pOdfi, Qj sin pOdfi , 

~i -i 

.and these may be made to depend on the values of the integrals 

\^ - jW 

-l -i 

It is proved that 

+ i 

Q£ sm K 0dfi - c ) ^ -; — \ — -f— -) -, — -f- , it n - <r be even, 

J (n - o-) !! (cr- A- 2) 11 (n + X + 1)11 

_i 

= if n - a- be odd, 

[iQ* sin A 6^/x - c <tt"l:^^^ if-ra-o-beodd, 

(n-o-- 1)11 (or- A- 2)11 (n+ A + 2)11 
-l 

= if n - o- be even. 

The factor c is equal to 2 or to it according as or + A is even or odd. 

5. The integrals 

+ i 



Q*Q? ( //x 



-1 

1 <a 
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are obtained in the form of a series having a finite number of 
terms. 

6. To find 

+i +i 



[q« 



Q*8mpdd[A and \ Q* cos pOdp, 
-i -i 

we may either express Qj or the trigonometrical functions in terms of 
a series of powers of sin 6. The second alternative leads to results 
which in general are more convenient. 
If we put 



G 



= 1; Oi = p; C 2 = p.£; C 3 =p.^L~l±- ; 

o = i^ + x r 2 iL ! - 

B =lj Bi=2>; B 2 =^-— -| ; B 3 = ^— z-^-^----- 

we find if <r be even, p odd, and w even, 

[r\ - aj (n + (T-l)\\(n~-3)\\ f n -, . -, 

Q sin pvdp = 7T \- ' ._A ^L__ <^ Ci . o- - 1 . o- + 1 

J •» ^ (?i - o-) ! ! (w + 2) ! ! L 

^ Or- 3.CT- 1 .0"+ 1 . O" +3 n (T — 5 . O" - 3 . O" — 1 . O" + 1 . (7 + 3 . O" + 5 

w - 3 . ?i + 4 ° n - 3 . n - 5 . n + 4 . n + 6 



« * » o 



and if cr be even, p even, and n even, 



TV 



(n-(r-l)!H»+3)!! L w-4.t? + 5 

-o o- — 5.a-—3.cr— l.cr+l.cr + 3.<r+5_ 

w - 4 . w - 6 . n + 5 . n + 7 
If n + o- + j? be even, the integral is zero. 

Similar equations are obtained for 

+ 1 
Q£ cospOdjji. 

-1 

7, The final results are expressed as follows : — 

If g% denote the coefficient of R^cos cr<f> in the series of spherical 
harmonics when E " = tfjQj, and ^ is a numerical coefficient, </> being 
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the longitude, the results of the investigation may be put into the 



form 



p — n + l 

g* = S a£ r * ? when n is even, o- odd, and p odd. 

jp = " 
p = n+l 

S << „ ti „ odd, o- „ „ j? even. 

•p = 91+1 

= 2 ^p^» j j w- ,, even, o- even, „ p odd. 



2> = 1 



■p = n+l 

= 2 &J s£ „ % „ odd, o- „ ., |? even. 

// — **> 

In these equations the factors a% b% are the coefficients of the 
Fourier series (see § 2), and the quantities r£, <", mj, si are 
numerical quantities, which (as well as their logarithms) are given in 
tables at the end of the paper as far as n = 12, <r = 12, p = 12. By- 
means of these tables the numerical work is reduced to a minimum, and 
the coefficients of the series may be obtained as far as terms of the 
12th degree. 

8. The proposed method is specially adapted to deal with problems 
like that of terrestrial magnetism, in which the function to be 
obtained as a series of spherical harmonics is not given directly, 
but by means of its differential coefficients. The force directed to 
the geographical north may by Fourier's analysis be obtained jas 
a sum, the terms of which have the form cos o-<£ cos pO, and sin o-^> 
cos jo0 when <t is even, and the form cos vcf> sin pO, sin cr<f> sin p6 when a- 
is odd. Integrating with respect to (9, the magnetic potential is obtained 
in a form such that the transformation into the series of spherical 
harmonics may be proceeded with. A separate expression of the 
magnetic potential is derived from the force directed to the geographical 
east. 
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